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ABSTRACT 

We generalize the effective potential to scalar field configurations which 
are proportional to the Hubble parameter of a homogeneous and isotropic 
background geometry. This may be useful in situations for which curvature 
effects are significant. We evaluate the one loop contribution to the Hubble 
Effective Potential for a massless scalar with arbitrary conformal and quartic 
couplings, on a background for which the deceleration parameter is constant. 
Among other things, we find that inflationary particle production leads to 
symmetry restoration at late times. 
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1 Introduction 



Consider two Heisenberg picture states, |\Efi) and \'^2), and a scalar field 
operator f{x) whose action is S[(p]. We define the background field ^{x) as 
the matrix element of (p{x) between these states, 



^1 



(f{x) 



(1) 



We define the quantum field as the difference between ip{x) and 

V9(x) = + . (2) 

The "effective field equation" for is obtained by taking the matrix ele- 
ment of the classical field equation. 
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6(p{x) 
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These definitions might seem a little pointless because the Heisenberg 
field ip{x) must evolve so as to make 5S/5(p{x) vanish. In that case ([3]) also 
vanishes and one would compute by substituting the solution for ^p{x) 
into (JT]). However, for many purposes it is simpler to compute perturbative 
operator corrections to for arbitrary then substitute these into 

([3]) and evaluate the matrix element assuming the matrix element for 
vanishes. The result is a C-number functional of $(a;) whose vanishing de- 
termines the actual background field. Proceeding in this manner is how one 
obtains the classic loop expansion for effective action r[$] [U 121 E], 



r[$] =5[$] + ^ln|det 
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+ 0{h^) . (4) 
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Here the exponent 



denotes the functional inverse, 



n -1 



/ 
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The scheme described above obviously depends a little upon the states 
For the familiar, in-out effective field equations the state |\l/2) is free 
vacuum (centered around $) in the infinite past, while is free vacuum 
(again centered on $) in the infinite future. For the Schwinger-Keldysh 
effective field equations one has = {"^2), and they are taken to be free 
vacuum (centered on $) in the infinite past (HE]. One can consider variants 
of the Schwinger-Keldysh equations in which the two states are still the same 
but the two states are defined at some finite time [HI E] . This is especially 
convenient for cosmological settings in which the initial singularity precludes 
reaching the infinite past, and we do not know what is a reasonable state to 
assume for the infinite future. 

The hard part about computing the effective action to a fixed loop order 
is solving for the functional inverses . Multiplying by i (times h which we 
henceforth set to one) gives the propagator in the presence of background $, 



An elegant formalism exists for deriving asymptotic expansions which are 
useful near coincidence, i.e., for x'^ near [51 El IS]- This allows one to give 
a wonderfully general treatment of ultraviolet divergences [TOl [HI [12] , but it 
sacrifices the nonlocal, ultraviolet finite contributions which constitute some 
of the most potentially interesting quantum effects. 

Solving for this propagator is typically so difficult that one is forced to 
specialize the effective action to a handful of special backgrounds. One of 
these special backgrounds, and one that plays an important role in flat space 
physics in trivial backgrounds, is $(x) = $0? where $0 is a constant in space 
and time. In that case the dependence of the kinetic operator on $0 typically 
reduces to a mass term and we can compute the propagator. 

The background $(x) = $0 defines the effective potential, 
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iA[$](x;x') . 
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In some cases V^ff($o) is precisely what we want. For example, if there is a 
vacuum state in flat space quantum field theory then it must be invariant 
under spacetime translations, so the actual solution is constant and we are 
guaranteed to find it by minimizing the effective potential. The value of the 
effective potential at this solution gives the vacuum energy density, and the 
second derivative of V^s at the solution determines the scalar mass. 

However, the more usual situation is that we would like to have the full 
effective action, and we only settle for the effective potential because it is 
what can be computed. There is often no point in working out the effective 
action for a constant background field versus some other configuration with 
spacetime dependence. It is especially pointless to demand that the back- 
ground be constant in cosmological settings, for which the geometry and 
many of the fields evolve in time. 

In this paper we shall explore a different class of backgrounds, namely 
those for which $(a;) is proportional to the Hubble parameter H[t), 



Of course the classical part of this Hubble Effective Potential — let us call 
it ^ub('^*o) — is straightforward to read off from the classical action. The 
derivative of the quantum correction to Vjiub is defined by the difference 
between the full effective field equations and the classical ones at the back- 
ground ([8]), 



As is apparent from the possibility of a conformal coupling in the classical 
Lagrangian, the dependence of V\mb{^oH) on H{t) is not limited to its ar- 
gument ^qH. We will see that quantum corrections to Vhub can also depend 
explicitly upon time. 

Knowing the Hubble Effective Potential might be useful in situations for 
which the scalar tends to evolve in tandem with the expansion of the universe. 
To see that this can occur in the early universe it suffices to study the classical 
equation of motion for a spatially homogeneous, quarticly coupled scalar on 
a homogeneous and isotropic background. 



= $0 X Hit) . 



(8) 





+ ?>H{t) + + iRit) + - = , 
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where m is the scalar mass, ^ the conformal couphng, A is the quartic couphng 
constant and R = 6H + 12H^ is the Ricci scalar. In the limit of large H the 
mass term becomes irrelevant. Let us also assume that the deceleration is 
constant, which is nearly the case during long epochs of cosmo logical history. 
In that case we have, 

H = -eH^ and H = 2e^H^ , (11) 

for some constant e. Under these assumptions — which certainly pertain at 
early enough times — substituting the Hubble scaling ansatz (jS]) into ffTUl) 
gives the following cubic equation for the proportionality factor $0) 

2e^ - 3e + 6(2 - e)^ 

It is not difficult to imagine values of ^ and e (e.g., power law inflation) 
for which the two symmetry-breaking solutions are real and the field tracks 
the evolution of the Hubble parameter. At least that would be the classical 
result; the Hubble Effective Potential describes how quantum corrections 
might change this. 

The usual Effective Potential Veff($o) tells us the energy density of the 
minimum uncertainty state centered on $o- What the Hubble Effective Po- 
tential gives us is that part of the energy density of the minimum uncertainty 
state centered on ^oH{t) which is not attributable to the kinetic energy of 
the scalar background. Of course this includes the usual result from Ves but 
it also receives contributions from the way the expansion of spacetime excites 
the scalar. When spacetime expansion leads to physical particle production 
there can be secular effects. 

This paper consists of five sections, of which the first is nearing its end. 
In section 2 we lay out the Feynman rules for using dimensional regular- 
ization to compute in massless Ay^^ theory on a homogeneous and isotropic 
background geometry of constant deceleration. These rules are employed in 
section 3 to obtain a fully renormalized result for the one loop correction to 
the Hubble Effective Potential. One fascinating feature of the result is sec- 
ular dependence in addition to that inherited from the simple fact that the 
scalar background is being evaluated at a time-dependent value ([8]). This is 
discussed in section 4, as is the dependence upon the arbitrary (but constant) 
value of the deceleration parameter. Our conclusions comprise section 5. 



:$n = . 



(12) 
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2 Feynman Rules 



The purpose of this section is to give the Feynman rules for computing loop 
corrections for a massless scalar with a quartic self-interaction on a nondy- 
namical background geometry with constant deceleration. We begin with a 
discussion of the geometry, then present the bare Lagrangian and the coun- 
terterms. The section closes with a discussion of the propagator in a general 
Hubble-scaled background ([8]). This is the key result which has made the 
present work possible. 

On the largest scales, the visible universe is well described by a homoge- 
neous, isotropic and spatially flat geometry. The invariant element takes a 
familiar form when expressed either in co-moving [t, x) or conformal [rj, x) 
coordinates, 

ds'^ = g^udx'^dx'^ = —dt'^ + a^{t)dx-dx = —drf + dx-dx . (13) 

To facilitate dimensional regularization we shall work in D spacetime dimen- 
sions, so the indices /i and v run from to {D — 1). 

Two derivatives of the scale factor a{t) have great importance, the Hubble 
parameter H{t) and the deceleration parameter 

, = -1-^^-1 +e(t). (14) 

Although the Hubble parameter may have changed by as much as 57 orders 
of magnitude from primordial inflation to now, the deceleration parameter is 
only thought to have varied from nearly —1, during primordial inflation, to 
+1, during the epoch of radiation domination. Because q{t) has been approx- 
imately constant for vast periods of cosmological evolution it is interesting 
to consider physics during an epoch of constant deceleration. 

It is simple to reconstruct the Hubble parameter and the scale factor 
over any epoch during which e is constant. If we define the initial values as 
i7(0) = Hq and a(0) = 1 then relation (fHll imphes, 

^W = r7^ a(t) = [l + ei/ot]' . (15) 

The conformal time rj is defined by drj = dt/a{t) up to a constant of integra- 
tion. We can choose this constant to give equivalent expressions for H and 



5 



H = — and a = — . (16) 

[-{l-e)Hov]— [-il-e)Hov] — 

Note that the universe is accelerating for < e < 1, and our conventions 
cause rj to approach zero fom below at late times. For e > 1 the universe is 
decelerating and r] is positive. Note also the simple and very useful relation, 



'l-e)Ha = — 



V 



(17) 



Finally, it is worth noting that the restriction to constant e reduces all cur- 
vature invariants to algebraic functions of e times powers of H, 



R 
UR 

R^^R^iu 



{D-l){D-2e)H^ , 
2{D-l)e{D-2e){D-l-3e)H^ , 

(D-l) \{D-1)D - A{D-l)e + De^ 



R'^'^'^Rpa,. = 2(D-1) 



D-Ae + 2e' 



The bare Lagrangian of interest for us is 
C = 



1 11 
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(18) 
(19) 

(20) 
(21) 



(22) 



Here y^o is the bare field and and Aq are the bare conformal and quartic 
coupling constants. We assume zero bare mass and, because mass is multi- 
plicatively renormalized in dimensional regularization, we shall not require a 
mass counterterm. 

The renormalized field ip is defined in the usual way. 



ip = Z 2 V9o • 

This gives the following expression for the bare Lagrangian ( 122|) . 



-9 



(23) 



(24) 



Renormalization is implemented by defining the bare couplings in terms of 
the physical ones plus counter parameters. 



Z = l + 5Z 



and 



Z^Xo = X + 6X. (25) 
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One then expresses the Lagrangian in terms of primitive interactions involv- 
ing the renormahzed fields and couplings, plus a series of counterterms. 



1 „„„ „ „„ , 1 o„ , 1 



d^y/^g^'^d^ - ^i?v^ - -A^^v/^ iA{x; x') = i6^{x - x') . (28) 



^6Zd,^d,^g^''^ - -5i^^R^ - -5\^^^ . (26) 

We consider ^ to be a free parameter of order one in A. The various counter 
parameters are of order, 

5Z = 0{X^) , 6^ = 0(A) and 6X = 0{X^) . (27) 

Therefore 6Z is not required at one loop while both 6^ and 6X are. 

Of course the formalism we have just reviewed is valid for any background 
geometry and any scalar field. If only one could do the same for the prop- 
agator! The equation it obeys is simple enough to write down for arbitrary 
metric and scalar backgrounds, 

1 

2' 

Alas, the set of known solutions is restricted to a pitiful handful of back- 
grounds. 

The recent technical advance which enables the present work is a tractable 
resolution for the problem of infrared divergences in the closely related propa- 
gator for the massless, minimally coupled scalar on a background of arbitrary, 
constant deceleration [13]. To make the connection, let us use conformal 
coordinates, and some facts about our constant deceleration geometry, to 
re-express the massless, minimally coupled kinetic operator, 

d^V^g^'-'d, = v^'Oy-'d,, (29) 

= a^~'{r^^^d,d^ + \[^)Ra'}a^~\ (30) 

= aT-^{r^^^d,d.+ ^ 4(l-e)V J ' ^^^^ 

= aT-i|r/'^-9^9, + ^^^H^}af-^ (32) 



Here the index Ummr is, 
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Multiplying by the factor {a' /a)~^^ — where a is the scale factor at and a' 
is the scale factor at x'^ — therefore carries the equation for the propagator 
of a massless, minimally coupled scalar, 

d^V^g^'d, zA^^,(a;; x') = , (34) 

to the form, 



z/2 -i 

mmc 4 



^^.^^^^^^m^J|(aa')T-i^A^^c(x;x')| =«r(x-x') . (35) 

A point of great importance for us is that equation fl35l) . with the re- 
placement z/jiimc ^1 can be solved for any constant index z/, whether or 
not z/ happens to obey ([33D [SI [13 [IS]- Now note that, for a scalar back- 
ground $(t, x) = ^oH{t), in a geometry of constant deceleration, the kinetic 
operator in (1281) can be reduced to precisely this form, 

= r/^"9^a^-29, - ^{D~l){D-2e)H^a'' - ^A^^/Z'a^ , (36) 
a^-^l'n^^-d.d, + , (37) 



7]- 



where the index v obeys. 



2 _ 1 /Z}-l-e\2 ^{D-l){D-2e) + iA<l>2 



^'--A—) (1-^)^ ■ '''' 

Hence we need only adapt the index of the previous solution [HI [151 [IB] once 
the infrared problem is solved. 

A formal solution to equation ( |35ll was found many years ago by Bunch 
and Davies for the spacetime manifold i?^ [H]. The generalization of this 
infinite space solution to D dimensions is [15], |16] , 



(47r) 2 i 



2 



^ D-l D w\ , , 
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where H and H' are the Hubble parameters at and x'^, respectively, and 



As stated, the solution appropriate for the massless, minimally coupled scalar 
has the index u = Urnmc, but solves our propagator equation (1251) as well, 
with the index u set to fl38|) . 

The infinite space solution fl39|) is only formal because it possesses infrared 
divergences in D = 4 [17] which the D-dimensional version inherits [18]. The 
presence of an infrared divergence in quantum field theory signals that the 
quantity which diverges is unphysical in some way. In this case the problem 
is that one is assuming the initial state could have been prepared in coherent 
Bunch-Davies vacuum out to infinite wavelengths. No causal process can 
accomplish this, so it is neither surprising nor problematic that making the 
assumption results in an infrared divergence. The assumption can be relaxed 
in either of two ways: 

1. One might continue to assume that the spatial sections are R^~^ but 
that the initially super-horizon wavelengths are in some less singular 
state [19]. That will guarantee there are no initial infrared singularities, 
which suffices to show that none develop later in time [2U] . 

2. One might instead work on the spatial manifold T^~^, which means 
there are no super-horizon wavelengths past a certain value and hence 
no possibility for an infrared divergence [21] . 

In [13] we obtained an explicit result for the latter fix. This makes the 
mode sum discrete, however, one can make the integral approximation with 
the lower limit of the co- moving momentum cutoff at some fixed value /cq. 
Then the corrections to the infinite space result fl39l) consist of the same 
integrand, integrated from = to \\k\\ = k^. Most of these corrections 
actually vanish for /cq = 0, and it is only necessary to include the ones which 
grow. To find them, note from (15^ that the infinite space result diverges 
whenever the index v in equation fl38|) obeys the relation. 





(40) 




(41) 



""^Note the unfortunate similarity between the geometrical parameter, e = 
the infinitesimal parameter e which enforces Feynman boundary conditions. 



H/H^, and 



9 



for some nonnegative integer A^. These correspond to the special values of 
V for which there is a logarithmic infrared divergence that dimensional regu- 
larization registers. The lower limit term which subtracts off this divergence 
is UJ, 



(47r)f r(i + z/)r(^) v-v^ 

x(A:oV) ' "'^CEE"'^^ll^^-^'ll'VV^^-'=-^ (42) 
with the coefficients ake defined as, 

^ r(^)r^(i-.) 

k\e\{N~k-e)\T{k + ^)T{i+l-u)T{N-k~i+l-i^)' ^ ' 

For > the infinite space mode sum also contains power law infrared 
divergences that dimensional regularization incorrectly removes by automatic 
subtraction [13J. These are removed by lower limit corrections which corre- 
spond to 5zAM(a^;a:') for < M < A^. All the other corrections vanish for 
/co = 0, so the full propagator for v near is, 

iA{x; x') ^ iAoo{x; x') + ^ SiAM{x; x') . (44) 

M=0 

We have shown [13] that this form agrees with results previously obtained 
for the special cases of e = for general D [22], and e = | for D = 4 |23]. 



3 Hubble Effective Potential at One Loop 

The purpose of this section is to compute and renormalize the one loop 
correction to effective field equation of massless (p^ theory for an arbitrary 
geometry of constant deceleration in the Hubble-scaled background ([8]), 

^t,x) = %x H{t) . (45) 



We begin by writing the one loop effective field equations in general, then 
specialize to constant deceleration and the Hubble-scaled scalar background. 
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This permits us to evaluate the coincident scalar propagator using our results 

dSni) and dmn). 

The Schwinger-Keldysh effective field equation for is obtained by 

taking the expectation value of the Heisenberg operator equation for ip{x), 



A 



6X 
"6 



■ if 



-9 ■ (46) 



One breaks the full field (p{x) up into its (D-number expectation value 
and the quantum field 



V9(x) = <l>(x) + 

At one loop the result is, 

6S[if] 



n 



where = 



A , ... „ , SX 



(47) 



<l>(^]|0'(x)|^])v^-5^$i?v^ (^'V^ + oih') . (48) 
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The effective field equation (148!) is valid for any metric and scalar back- 
ground, but it is of course impossible to obtain explicit results for the coin- 
cident scalar propagator which is implicit in (f2|0^(x)|f2), except in special 
backgrounds. We therefore specialize at this point to the geometry of arbi- 
trary but constant deceleration described in the previous section, and to the 
Hubble-scaled scalar ([8]), 

{Q\<f)\x)\Q) = iA{x;x) + 0{ff) , (49) 

N 

= iA^{x; x) + 6iAM{x; x) + 0{h^) . (50) 

M=0 

Of the various one loop contributions to (l50l) . only the infinite space 
propagator harbors ultraviolet divergences, so we will evaluate it, renormalize 
and then include the finite contributions from 5iAM{x:,x). The coincidence 
limit of the infinite space propagator (15^ is. 



iAoo(a;;x) 



■^)r(i-f) 



(4vr) 



r(i+z/)r(i-z/) 



(51) 
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At this point we define S = A — D and expand the i/-dependent ratios, 



D-l 



v{\+un\-u) 



r 



D-3 



D-3 



(52) 



+ + (53) 



x<^ 1 



The expansion of the divergent Gamma function is also important, 



---1+7 + 0(5) 



(54) 



where 7 ~ .577216 is the Euler constant. Substituting ( !53ll54l) in (ISTI) gives. 



r(i-f) 



(47r)T[(l-e)2i72] 
1 



^6/ 2 . 



+ 



167r2 



+ 



.^•(i+-) + <i- 
-^i(l-5e+3ai/2 + OW. (66) 



Employing (!50|) and (!55|) in the effective field equation (148!) shows that 
we must make the following choices for the one loop counter parameters. 



5^ 
6X 



QJ (4 



2 (47r)2/i^ 



(56) 
(57) 



Here fii and /i2 are two dimensional regularization mass scales which can 
be used to subsume the finite parts of the counter parameters. With these 
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choices we can take the unregulated hmit to obtain Vl^^{^), 





A A ^ 

+7:7r^ {l-5e + 3e^)H^^ + o 5Z ^^^m{x; x)^ + 0{ff) . (58) 



967r2 



Af=0 



A result for this has already been reported in the literature for the de Sitter 



case of e = and z/ 



1 /3-e\2 



Note that for D = A the index fl38ll becomes 



D=4 4 V 1 - 



[^(12-66) + lA$2] 1 [{^-\)R+\W^ 



1 



(59) 



We will assume v is positive (corresponding to inflation with small $o) but 
one should bear in mind that it might be negative or even imaginary. 

It remains to evaluate the lower limit contributions 5iAM{x\x) which 
were defined in equations (14211431) . At coincidence we have ||x — = and 
7] = rj'. Only the k = term of the sum in expression fH2l) contributes, and 
the sum over i can be performed to give, 



6iAM(x: x) 



327r2 



r(2z/-M)r(2i/-2M) 

r(M+i)r2(i/+i-M) ■ 



(60) 



It might be worth noting from relation (fT7|) that the factor k^rf can be 
written in terms of the Hubble parameter and scale factor. 



2„2 



'1- 



(61) 



Note also that if the sum over M is extended to infinity the result can be 
expressed hypergeometric function. 



m[x;x) 



M=0 



27r3 



2 



X 



T\u)2F^{^--u,'^--u-l-2u,l-u,^-u--kW) . (62) 
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When the index v approaches | + A^, for some nonnegative integer N it 
is desirable to extract the infrared divergence from ?/'(| — z/) in expression 



1 e; ^ 1 



(63) 



M=-\ 



A manifestly finite form can be obtained by combining the series in f l63|) with 
the series of lower limit contributions. This is facilitated by the identity, 



(64) 



The final result for is 
A 



Kub(<f ) = 

4. 



967r2 



l-5e + 3e2)i72$ 



+ < 4 + 



327r2 
A2 



In 



6 647r2 
l-e)2ij2 



In 



327r2 



N 



M=-l 2 



T{2iy-M)T{2u-2M) {Aklr] 
2r(M+l)r2(z/+i-M) 



(^'-^)|a$ + 0(;i^).(65) 



4 Field Strength, Time and Deceleration 

In the last section we obtained explicit forms (!58l) and (!65|) for V^'ub('^')- The 
purpose of this section is to discuss how the Hubble Effective Potential de- 
pends upon the scalar field strength, upon time, and upon the arbitrary con- 
stant deceleration. We begin by deriving a large field expansion and checking 
the flat space correspondence limit against the classic result of Coleman and 
Weinberg ^25j . This suggests symmetry breaking, which we confirm, for early 
times, by expanding around $ = 0. However, the explicit time dependence of 
the Hubble Effective Potential introduces a profound change: prolonged time 
evolution during inflation restores the $ — $ symmetry. We develop an 
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expansion for intermediate field strengths which shows that the symmetry 
breaking minima come before this point. Then we investigate the exphcit 
time dependence more fully, and the section closes with a discussion of the 
dependence upon e. In all cases we assume < e < 1, corresponding to 
inflation. We shall also equate the two mass scales /ii = /i2 = A* omit 
the +0{h'^) at the end of each result for ^hubl*^) order to obtain more 
compact expressions. 

When the scalar field strength is large the index u becomes imaginary, 



[(e-|)i?+|A<i.^ 



1- 



2^2 



1 

4 



(66) 



In this case it is best to expand the first form ( l58l) . The key relation we need 
comes from differentiating the logarithm of Stirling's formula for the Gamma 
function, 

V^(^) = H^) -l-lL + T^ + 0(^). (67) 



2z 12^2 I20z^ 
The resulting expansion of the digamma function terms is. 



In 



(e-i)/?+iA<|.^ 



(l-e)2/72 



(e-i)/?+|A$2 



1 

^ 15 



n 2 



(e-i)i?+|A$2 



0{^-'') . (68) 



When the index is imaginary there are no infrared divergences at all and 
the lower limit correction terms can be dropped. From (!58|) we see that the 
resulting large $ form is. 



A 

1 



A 



A 



967r2 
A 



^6/ 2 



l-5e + 3e^)H'^<l> 

(e-i)i?+iA$^ 



In 
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(e-i)i?+lA<|.2 



+ 0($-'') V<l> . (69) 



The coefficient of the logarithm term agrees, as it must, with equation (3.10) 
of the classic paper by Coleman and Weinberg |25j. It agrees as well with the 
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result previously obtained in de Sitter [21] and with what has been shown in 
general backgrounds . This large field strength regime is also appropriate 
to a Higgs with ^ S> 1, which has been proposed as a possible candidate for 
the inflaton [271128]. 

Expression fl69l) suggests that symmetry breaking may occur, however, 
care must be taken of the fact that fl5^ is only the large field form of a 
function that may not actually pass through zero for $ 7^ 0. Symmetry 
breaking in the large field form (l69l) requires the logarithm to be negative, 
which is realized for, 

(e-^)i?+^$^</^^ (70) 
On the other hand, the large field expansion is valid for, 

{^-l)R+^<^'>il-e?H'. (71) 

These two conditions can only be consistent if, 

(1-6)^/72 < ^2 _ ^72) 

This is very dubious if the theory is renormalized at some low energy scaled 
To check if symmetry breaking really occurs we need to expand ^ub(^) 
around $ = 0, and the best form for this is expression fl65|) . For small 
coupling (A <^ 1) we can see right away that the sign of the conformal 
coupling controls whether or not the potential dips below zero, 

vUi^) = ^m + oiX). (73) 

For ^i? < there will be symmetry breaking, whereas there will be no sym- 
metry breaking (at least near $ = 0) if ^i? > 0. The only case requiring 
detailed analysis is minimal coupling (^ = 0), or near minimal coupling. 

^ Of course one could use the curved space renormalization group 29J to move the scale 
up to the level of primordial inflation. We could not set /z to the instantaneous Hubble 
parameter H{t) because that is time dependent [30], but it would be perfectly possible to 
evolve /i up to the constant Hq. Then the large logarithms ln(7J^//i^) that appear in our 
results would be absorbed into redeflnitions of the conformal and quartic couplings, ^{fJ.) 
and A(/x). As long as the low scale A is small enough to dominate the large logarithms the 
net effect would be no change, so we shall continue to imagine that /x is a Standard Model 
scale very much smaller than the Hubble parameter. One important consequence is that, 
when we assume "minimal coupling," we mean f (/i) ~ at Standard Model scales. 
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When ^ = 0, the value of the index z/ at $ = is, 



= hm 1/ 



1 /3-e 



2 Vl-e 



(74) 



We shall assume < e < 1, which implies z/q > f- Let us also assume that 
uq lies in the range, 



iV+l<z/o<iV + 2, 



(75) 



for some nonnegative integer N. Whether or not the potential falls below 
zero near $ = is controlled by the second derivative. 



limW' (0) = ^ 



X • 



i2-e) 



In 



(5-46) 



N 



M=0 



^ Uq-M- 



2 1 



r(2i/o-M)r(2i/o-2M) 



° 4 2M!r2(z/o + |-M)(4fc2^2; 



• (76) 



The sign of expression (17^ is basically determined by the competition 
between its first and last terms, 



ImCb(O) 



AiJ2 
32^ 



-(2-e) In 



;i-e)2ij-2 



(l-e)2 8r2(z.o) / 4 



^0-2 



k^rf 



(77) 

The first (logarithm) term tends to make Vjjub(O) < 0, whereas the second 
('^o ~ f pole) term tends to make it positive. For a renormalization scale 
yU which is appropriate to the Standard Model, and a much higher scale 
Hubble parameter, the first term is large. For an infrared cutoff of /cq = -f^O) 
corresponding to no super-horizon modes in the initial state, the product 
k^rf is initially of order one, but rapidly approaches zero. This means that 
the second term is initially smaller in magnitude than the first, but grows 
without bound at late times. Hence we expect symmetry breaking shortly 
after the initial time of rj^ = —l/{l—e)HQ, with symmetry restoration at late 
times. Because the exponent z^o — | grows with e, we also expect the onset of 
symmetry restoration to come sooner for larger values of e. Both expectations 
are born out by explicit numerical evaluation of the full expression fl76|) — 
not just the two terms of fl77|l — as displayed in Fig. [H The symmetry 
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Figure 1: Both plots show ^^"^(0) ^ Aif^/327r^ as a function of e for ^ = 
with (i/Hq = 10~^^ and ko/HQ = 1. The left plot is for HqI] = —1 (shortly 
after the initial time) and the right plot is for the later time HqI] = — 10~^. 
Note that all values of e within the plotted range show symmetry breaking 
initially, and that time evolution tends to restore the symmetry. Symmetry 
restoration takes place sooner for larger e. 



\ 0.015 




\ 0.010 




\ 0.005 




-3 i--'^''^ 





Figure 2: Plot of the Vq = Vhub/-f^^ as a function of $o for ^ = 0, A = .01 with 
e = .1, fi/Ho = 10~^^ and ko/Ho = 1. The left plot is for HqI] = —1 (shortly 
after the initial time) whereas the right plot is for HqI] = —10"^, correspond- 
ing to about 12.8 e-foldings of further evolution. Symmetry restoration occurs 
at about HqI] ^ —1.257 x 10^^ for this value of e. 
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restoration effect can also be seen by evolving the full potential for a fixed 
value of e, as is displayed in Fig. [21 

We have seen how to obtain expansions for the large field and small field 
regimes. As the magnitude of $ increases, the value of v decreases to zero, 
and then becomes imaginary. It is possible to develop an expansion for the 
intermediate regime of $, in which v is close to |. The small quantity in this 
case is. 



[(e-|)i?+|A$^] 



(78) 



Regarding ( ITHI) as small amounts to expanding in powers of the parameter 
IS.V defined as, 



1 

Ai/ = V 

2 



(e-|)i?+|A$^ 

fl-e)2/72 



+ 



-\ 2 



:i-e)2iJ2 



(e-|)i?+|A$2 



-\ 3 



'l-e)2iJ2 



+ 



(79) 



The digamma functions have the following expansions, 

oo 

n=2 
^ oo 



(80) 
(81) 



n=2 



where is the Riemann zeta function. It is also useful to note. 



(82) 



Because v is less than the = pole at = |, we can ignore the infrared 
corrections (which fall off and are also suppressed by a prefactor proportional 
to Az/^) and the result is. 



1 2 A 

18^3^~5767r2 



1 1 



A(l-e)2i72Az/(l-Ai/) 
32^^2 



In 



i?+ (l-e)2i/2Az/|l-^Az/4 
l-e)2i72. 



A 



3 3' 327r2 

oo 

27-2^C(2n+l)A 



2n 



n=l 



(83) 
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For Az/ = the result has the same sign as $ for ^ = and A < 32-7?^, so we 
are beyond the symmetry breaking minima (if any), although still not in the 
large field regime. 

We turn now to a systematic investigation of explicit time dependence. 
This enters the Hubble Effective Potential in two ways: 

1. From the factors of ln[if (t)] that combine with the logarithms of the 
dimensional regularization mass scale /xi = = A* in expressions (155|) 
and fl^ : and 

2. From the factors of (kQ-r]'^)^^'^^''^ that multiply the infrared corrections 
of expression (!65l) . 

The first of these has been included in the large field expansion (!69|) and in 
the small Az/ expansion (!83|) . If the field $(x) plays a role in low energy 
physics then the scale fi is presumably in the Standard Model, and hence 
much less than the Hubble parameter H(t) at all times during primordial 
inflation. The Hubble parameter does not change much during inflation so 
we expect these logarithms to be large and slowly decreasing, which reduces 
the coefficients of the i?$ and $^ contributions to ^ub(^)- have seen 
from expressions ( !69|) and ( !83l) that this time dependence has little effect on 
symmetry breaking. 

The second source of temporal dependence behaves very differently. Re- 
call that the conformal time t] begins at rjo = — 1/(1 — e)ifo during inflation 
and approaches zero at asymptotically late times. This means that the fac- 
tors of (A;q?7^)^^"'"2~'^ grow like powers of the scale factor for any M such that 
M < — \- The fastest growth comes from M = 0, which contributes the 
following term to V^ubl*^); 

The initial factor of H^{t) actually makes (18^ fall with time, but it still grows 
with respect the H^^ terms that would otherwise drive symmetry breaking. 
Another point to note is the dependence on the constant $o- Because the 
index v decreases as a function of eventually reaching zero and becoming 
imaginary, we see that flMl) vanishes in the large fleld limit. The same is true 
for the contributions from any M. 

The various trends we have just described can be seen by plotting the 
Hubble Effective Potential for the same parameters at different times. To 
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factor out the time dependence of the Hubble parameter we have scaled the 
potential and the field by H(t), 



For the various parameters we chose, 

e = , A = .l , e=^ , ^i = 10-''Ho , ko = Ho. (86) 

For these values, the transition from symmetry breaking to symmetry restora- 
tion occurs at about HqI] ^ —5.05 x 10~^. Fig. [3] shows Vq over the range 
—2 < $0 < +2 for three different conformal times, 

HoV = -lO-' , HoV = -lO'' , HoV = -lO-'. (87) 

The earliest time shows symmetry breaking, whereas the symmetry has been 
restored by the second time, and the third plot shows that further time 
evolution makes the potential steeper at the origin. Fig. H] shows the same 
three curves but over the expanded scale —8 < $0 < +8. For parameters 
( !86ll the center of the intermediate regime — expansion (1831) — comes at 
$0 = i6, and the index u goes imaginary at $0 = ±14/^5 ^ ±6.3. The 
large field expansion fl69|) is valid by $0 = ±8. Note that the steepening at 
late times is a small field effect; time evolution has very little impact in the 
large field or even the intermediate field regime. 

We close with a comment on e dependence. This is only very weak in the 
large field regime (169|1 . as required by the need to agree with the result of 
Coleman and Weinberg [25]. On the other hand, e dependence is an order 
one effect in the intermediate regime (I83|) — recall that R = Q{2-e)H^ and 
that Au vanishes for $q = 2(2 — e)/ A. In the small field regime the value 
of e controls the rate at which time evolution causes the M = infrared 
contribution fl8^ dominate. One can see all this by plotting Vq versus $0 for 
different values of e at the same time. In Fig. [5] and Fig. [6] we have chosen, 

^ = , A = .l , /x = 10-^3ifo , ko = Ho , HoV = -lO-\ (88) 

for three different values of e, 

e = .1 , e = .3 , e = .4 . (89) 
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Figure 3: Three different snapsfiots of Vq = Vhnh/H'^ as a function of $o for 
^ = 0, A = .1 with e = .2, fi/Ho = 10'^^ and ko/Ho = 1. The red dashed 
curve is for HqI] = — 10~^, the purple dotted curve is for HqT] = —10"'^ (about 
3.7 e-foldings later), and the blue solid curve is for Horj = — 10~^ (another 
3.7 e-foldings later). 



22 




Figure 4: The same three snapshots of Vq = V\mb/H^ as a function of $o 
as in Fig. [3] over a larger range of $o- For these parameters (A = .1 and 
e = .2) the quantity Az/ = | — z/ of the intermediate expansion ( |83l) vanishes 
at $0 = i6. The index u becomes imaginary for $o = ~ ±6.3 and the 
large field expansion (1^ is valid by the end of the plotted range at $o = ±8. 
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Figure 5: Vq = VhubZ-f^"^ as a function of $o for three different values of e. 
The other parameters are ^ = 0, X = .1, fi/ Hq = 10~^^, ko/Ho = 1 and 
HqI] = —10^^. The red dashed curve is for e = .1, the purple dotted curve is 
for e = .3, and the blue solid curve is for e = .4. 

In Fig.[5]the three curves are plotted over the range —2 < $o < +2- The plot 
for e = .1 shows symmetry breaking, whereas the symmetry has been restored 
for e = .3, and the plot for e = .4 shows that the potential steepens with 
increasing e. In Fig. [6] the same three curves are plotted over the expanded 
range —8 < $o < +8- For these parameters the center of the intermediate 
regime (183|) is at about $o ~ ±6, and the large field regime (169|) pertains 
by $0 = =t8. One can see that the curve for e = .4 is still noticeably above 
the other curves in the intermediate regime, but the three curves are almost 
indistinguishable in the large field regime. 

5 Discussion 

We have introduced a new sort of effective potential which might be use- 
ful in cosmological settings for which the curvature dominates over other 
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Figure 6: The same plots of Vq = Vauh/H^ as a function of $o as in Fig. [5] 
over a larger range of $0- For these parameters the center of the intermediate 
regime (l83l) is at about $o ~ =t6, whereas the large field result (l69l) applies 
by $0 = ±8. 



25 



dimensionful parameters. Instead of evaluating the effective action for a con- 
stant background scalar, we evaluate it at a background that scales with the 
Hubble parameter, 

^{t,x) — >^oy<H{t). (90) 

In section 2 we show how to obtain the propagator (which is the limiting 
step on any computation of the effective action) for a massless scalar with 
arbitrary conformal coupling, in an FRW geometry of arbitrary but constant 
deceleration. In section 3 we evaluated the Hubble Effective Potential at one 
loop order for a quartic self-coupling. 

Section 4 was spent discussing the dependence upon field strength, time 
and deceleration under the assumption of primordial inflation far above Stan- 
dard Model scales. We obtained expansions valid for large field strength fl^ 
and for intermediate field strength fl55]) . and we made a careful study of the 
curvature (176!) at $o = 0. Of course the large field limit agrees with the 
fiat space result of Coleman and Weinberg [25], which suggests that sym- 
metry breaking may occur. We confirmed that it does occur, but only at 
early times. As time progresses the $ — $ symmetry is restored, with the 
effect coming sooner the larger the deceleration. Subsequent time evolution 
steepens the potential near the origin, without of course affecting its large 
field limit. 

The fascinating time dependence of the Hubble Effective Potential seems 
to mostly derive from inflationary particle production. (There is also a small 
effect due to the time dependent Hubble parameter declining relative to the 
fixed scale of renormalization.) One can understand what is happening on a 
simple, qualitative level by noting that inflationary particle production will, 
of necessity, increase the field strength. So it drives the field away from 
the small symmetry breaking minima. This is why the effective potential 
is driven up, at small field strength, but not affected much at large field 
strength. 

One consequence of the explicit time dependence we have found is that 
the effective field $(t) cannot really evolve according to the Hubble-scaling 
ansatz which made the calculation possible! This does not mean what 
we have done is irrelevant, any more than it would be irrelevant to use the 
effective potential to show that the field cannot be constant. We now have 
a proof thai the field cannot exactly track H{t). Further, as long as the 
time dependent terms are small, the Hubble Effective Potential probably 
describes the actual evolution well, just as the effective potential is assumed 
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to do for slowly varying fields. At the very least, one can use our formalism 
to determine which way the actual background is pushed away from exact 
Hubble-scaling. 

The effective potential we obtain is quite similar to those derived on de 
Sitter background for this same model [21], for Yukawa theory [31] and for 
scalar quantum electrodynamics (SQED) [32j. One striking feature of those 
models is that the effective potential which describes scalar evolution is not 
quite the same as the one that appears in the effective stress tensor and would 
govern the gravitational response if back-reaction were included. This can 
have interesting consequences. For example, although the effective potential 
for SQED is minimized at zero field strength, the vacuum energy actually 
decreases, for a while, as the field strength increases [32] • The reason for 
the differences was traced to the fact that the effective potential depends 
nontrivially upon the metric, so its variation — which gives the stress tensor 
— must reflect this dependence [311 [32] . The same thing should be true for 
our model, and it would be quite interesting to quantify this expectation by 
working out the one loop correction to the gravitational field equations. It 
would be particularly interesting to explore the gravitational consequences 
of symmetry restoration. 

Finally, we should contrast our view of the infrared issue with some differ- 
ent approaches that have been suggested to dealing with infrared divergences 
in loop corrections to the power spectrum [33] • One might argue that be- 
cause super-horizon modes have no apparent spatial variation, they are not 
observable, and the results of quantum field theory must be changed to ex- 
clude them. It seems to us that spatial constants are observable — indeed, 
the current universe seems to be experiencing the effect of something very 
much like a small cosmological constant [34j ! 

We also believe that the infrared divergences of loop corrections to the 
power spectrum are clear analogs of the old problem which afflicts the prop- 
agator of a massless, minimally coupled scalar on many cosmological back- 
grounds [17|. In that case the infrared divergence results, as infrared diver- 
gences always do in quantum field theory, from posing an unphysical question. 
The physically dubious feature of the question is that one can prepare co- 
herent Bunch-Davies vacuum over an infinite spatial section. Finite results 
are obtained when this assumption is relaxed, either by working on a com- 
pact space [21] or by making the super-horizon modes of the initial state less 
singular than Bunch-Davies [19]. This sort of resolution leads to finite re- 
sults which, however, grow with time as more and more initially sub-horizon 
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modes are redshifted to long wave lengths [HI dSl [221 ES]- 

The fascinating secular effect of symmetry restoration we have found de- 
rives from precisely this source, and it seems to have observable consequences. 
For example, suppose we consider the scalar to be the Standard Model Higgs, 
and suppose that it is effectively massless and minimally coupled during pri- 
mordial inflation. During the early stages of inflation, when the symmetry 
is broken, the fermions of the Standard Model will be massive, with the ex- 
pected mass of the order ~ yH^o, where y is the Yukawa coupling and 
$0 ~ 1- When < H (which is valid for most Standard Model fermions 
with y <^1) they will be produced during inflation, reaching the occupation 
number [55] . 

7i , ~ \ \ ~ 1/2 f9ll 

^ exp(27rm^/i7) + 1 exp(27r?/$o) + 1 ' ' ^ ' 

This is a good estimate as long as H changes adiabatically slowly in time, 
i.e. when e -C 1. (For a finite e this formula will probably be corrected to, 

= l/[exp(27rm^/[(l — e)H]) + 1].) Now, if inflation is of a relatively short 
duration, the symmetry may be restored quite late during inflation, in which 
case some of the fermions may survive and carry an observable information 
about the duration of inflation. 

Another source of potentially observable effects is in cosmological pertur- 
bations. The effective potential we have been computing presumably affects 
gravity — although probably not quite the way a classical scalar potential 
would [311 [32]. That will change the geometry and alter the spectrum of 
density perturbations. And note that, whereas quantum gravitational cor- 
rections are necessarily suppressed by a factor of the small loop counting 
parameter GH^ < 10~^^, the effects of quantum matter loops need only be 
suppressed by Standard Model couplings. 
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